Range-wide and temporal genomic analyses reveal the consequences of near-extinction in Swedish moose.

Ungulate species have experienced severe declines over the past centuries through overharvesting and habitat loss. Even if many game species have recovered thanks to strict hunting regulation, the genome-wide impacts of overharvesting are still unclear. Here, we examine the temporal and geographical differences in genome-wide diversity in moose (Alces alces) over its whole range in Sweden by sequencing 87 modern and historical genomes. We found limited impact of the 1900s near-extinction event but local variation in inbreeding and load in modern populations, as well as suggestion of a risk of future reduction in genetic diversity and gene flow. Furthermore, we found candidate genes for local adaptation, and rapid temporal allele frequency shifts involving coding genes since the 1980s, possibly due to selective harvesting. Our results highlight that genomic changes potentially impacting fitness can occur over short time scales and underline the need to track both deleterious and selectively advantageous genomic variation.

Tradeoffs between resources and risks shape the responses of a large carnivore to human disturbance.

Wide-ranging carnivores experience tradeoffs between dynamic resource availabilities and heterogeneous risks from humans, with consequences for their ecological function and conservation outcomes. Yet, research investigating these tradeoffs across large carnivore distributions is rare. We assessed how resource availability and anthropogenic risks influence the strength of lion (Panthera leo) responses to disturbance using data from 31 sites across lions' contemporary range. Lions avoided human disturbance at over two-thirds of sites, though their responses varied depending on site-level characteristics. Lions were more likely to exploit human-dominated landscapes where resources were limited, indicating that resource limitation can outweigh anthropogenic risks and might exacerbate human-carnivore conflict. Lions also avoided human impacts by increasing their nocturnal activity more often at sites with higher production of cattle. The combined effects of expanding human impacts and environmental change threaten to simultaneously downgrade the ecological function of carnivores and intensify human-carnivore conflicts, escalating extinction risks for many species.

Mountain- and brown hare genetic polymorphisms to survey local adaptations and conservation status of the heath hare (Lepus timidus sylvaticus, Nilsson 1831).

We provide the first whole genome sequences from three specimens of the mountain hare subspecies the heath hare (Lepus timidus sylvaticus), along with samples from two mountain hares (Lepus timidus timidus) and two brown hares (Lepus europaeus) from Sweden. The heath hare has a unique grey winter pelage as compared to other mountain hares (white) and brown hares (mostly brown), and face regional extinction, likely due to competitive exclusion from the non-native brown hare. Whole genome resequencing from the seven hare specimens were mapped to the Lepus timidus pseudoreference genome and used for detection of 11,363,883 polymorphic nucleotide positions. The data presented here could be useful for addressing local adaptations and conservation status of mountain hares and brown hares in Sweden, including unique subspecies.

Large mammal telomere length variation across ecoregions.

BACKGROUND: Telomere length provides a physiological proxy for accumulated stress in animals. While there is a growing consensus over how telomere dynamics and their patterns are linked to life history variation and individual experience, knowledge on the impact of exposure to different stressors at a large spatial scale on telomere length is still lacking. How exposure to different stressors at a regional scale interacts with individual differences in life history is also poorly understood. To better understand large-scale regional influences, we investigated telomere length variation in moose (Alces alces) distributed across three ecoregions. We analyzed 153 samples of 106 moose representing moose of both sexes and range of ages to measure relative telomere lengths (RTL) in white blood cells. RESULTS: We found that average RTL was significantly shorter in a northern (montane) and southern (sarmatic) ecoregion where moose experience chronic stress related to severe summer and winter temperatures as well as high anthropogenic land-use compared to the boreal region. Our study suggests that animals in the northern boreal forests, with relatively homogenous land use, are less disturbed by environmental and anthropogenic stressors. In contrast, animals in areas experiencing a higher rate of anthropogenic and environmental change experience increased stress. CONCLUSION: Although animals can often adapt to predictable stressors, our data suggest that some environmental conditions, even though predictable and ubiquitous, can generate population level differences of long-term stress. By measuring RTL in moose for the first time, we provide valuable insights towards our current understanding of telomere biology in free-ranging wildlife in human-modified ecosystems.

A century of genetic homogenization in Baltic salmon-evidence from archival DNA.

Intra-species genetic homogenization arising from anthropogenic impacts is a major threat to biodiversity. However, few taxa have sufficient historical material to systematically quantify long-term genetic changes. Using archival DNA collected over approximately 100 years, we assessed spatio-temporal genetic change in Atlantic salmon populations across the Baltic Sea, an area heavily impacted by hydropower exploitation and associated with large-scale mitigation stocking. Analysis was carried out by screening 82 SNPs in 1680 individuals from 13 Swedish rivers. We found an overall decrease in genetic divergence and diminished isolation by distance among populations, strongly indicating genetic homogenization over the past century. We further observed an increase in genetic diversity within populations consistent with increased gene flow. The temporal genetic change was lower in larger wild populations than in smaller wild and hatchery-reared ones, indicating that larger populations have been able to support a high number of native spawners in relation to immigrants. Our results demonstrate that stocking practices of salmon in the Baltic Sea have led to the homogenization of populations over the last century, potentially compromising their ability to adapt to environmental change. Stocking of reared fish is common worldwide, and our study is a cautionary example of the potentially long-term negative effects of such activities.

African and Asian leopards are highly differentiated at the genomic level.

Leopards are the only big cats still widely distributed across the continents of Africa and Asia. They occur in a wide range of habitats and are often found in close proximity to humans. But despite their ubiquity, leopard phylogeography and population history have not yet been studied with genomic tools. Here, we present population-genomic data from 26 modern and historical samples encompassing the vast geographical distribution of this species. We find that Asian leopards are broadly monophyletic with respect to African leopards across almost their entire nuclear genomes. This profound genetic pattern persists despite the animals' high potential mobility, and despite evidence of transfer of African alleles into Middle Eastern and Central Asian leopard populations within the last 100,000 years. Our results further suggest that Asian leopards originated from a single out-of-Africa dispersal event 500-600 thousand years ago and are characterized by higher population structuring, stronger isolation by distance, and lower heterozygosity than African leopards. Taxonomic categories do not take into account the variability in depth of divergence among subspecies. The deep divergence between the African subspecies and Asian populations contrasts with the much shallower divergence among putative Asian subspecies. Reconciling genomic variation and taxonomy is likely to be a growing challenge in the genomics era.

Large-scale spatial variation of chronic stress signals in moose.

The physiological effects of short-term stress responses typically lead to increased individual survival as it prepares the body for fight or flight through catabolic reactions in the body. These physiological effects trade off against growth, immunocompetence, reproduction, and even long-term survival. Chronic stress may thus reduce individual and population performance, with direct implications for the management and conservation of wildlife populations. Yet, relatively little is known about how chronic stress levels vary across wild populations and factors contributing to increased chronic stress levels. One method to measure long-term stress in mammals is to quantify slowly incorporated stress hormone (cortisol) in hair, which most likely reflect a long-term average of the stress responses. In this study, we sampled 237 harvested moose Alces alces across Sweden to determine the relative effect of landscape variables and disturbances on moose hair cortisol levels. We used linear model combinations and Akaike's Information Criterion (corrected for small sample sizes), and included variables related to human disturbance, ungulate competition, large carnivore density, and ambient temperature to estimate the covariates that best explained the variance in stress levels in moose. The most important variables explaining the variation in hair cortisol levels in moose were the long-term average temperature sum in the area moose lived and the distance to occupied wolf territory; higher hair cortisol levels were detected where temperatures were higher and closer to occupied wolf territories, respectively.

Droplet digital PCR assays for the quantification of brown trout (Salmo trutta) and Arctic char (Salvelinus alpinus) from environmental DNA collected in the water of mountain lakes.

Classical methods for estimating the abundance of fish populations are often both expensive, time-consuming and destructive. Analyses of the environmental DNA (eDNA) present in water samples could alleviate such constraints. Here, we developed protocols to detect and quantify brown trout (Salmo trutta) and Arctic char (Salvelinus alpinus) populations by applying the droplet digital PCR (ddPCR) method to eDNA molecules extracted from water samples collected in 28 Swedish mountain lakes. Overall, contemporary fish CPUE (catch per unit effort) estimates from standardized survey gill nettings were not correlated to eDNA concentrations for either of the species. In addition, the measured environmental variables (e.g. dissolved organic carbon concentrations, temperature, and pH) appear to not influence water eDNA concentrations of the studied fish species. Detection probabilities via eDNA analysis showed moderate success (less than 70% for both species) while the presence of eDNA from Arctic char (in six lakes) and brown trout (in one lake) was also indicated in lakes where the species were not detected with the gillnetting method. Such findings highlight the limits of one or both methods to reliably detect fish species presence in natural systems. Additional analysis showed that the filtration of water samples through 1.2 µm glass fiber filters and 0.45 µm mixed cellulose ester filters was more efficient in recovering DNA than using 0.22 µm enclosed polyethersulfone filters, probably due to differential efficiencies of DNA extraction. Altogether, this work showed the potentials and limits of the approach for the detection and the quantification of fish abundance in natural systems while providing new insights in the application of the ddPCR method applied to environmental DNA.

Carnivores, competition and genetic connectivity in the Anthropocene.

Current extinction rates are comparable to five prior mass extinctions in the earth's history, and are strongly affected by human activities that have modified more than half of the earth's terrestrial surface. Increasing human activity restricts animal movements and isolates formerly connected populations, a particular concern for the conservation of large carnivores, but no prior research has used high throughput sequencing in a standardized manner to examine genetic connectivity for multiple species of large carnivores and multiple ecosystems. Here, we used RAD SNP genotypes to test for differences in connectivity between multiple ecosystems for African wild dogs (Lycaon pictus) and lions (Panthera leo), and to test correlations between genetic distance, geographic distance and landscape resistance due to human activity. We found weaker connectivity, a stronger correlation between genetic distance and geographic distance, and a stronger correlation between genetic distance and landscape resistance for lions than for wild dogs, and propose a new hypothesis that adaptations to interspecific competition may help to explain differences in vulnerability to isolation by humans.

Discovery of SNPs for individual identification by reduced representation sequencing of moose (Alces alces).

Monitoring of wild animal populations is challenging, yet reliable information about population processes is important for both management and conservation efforts. Access to molecular markers, such as SNPs, enables population monitoring through genotyping of various DNA sources. We have developed 96 high quality SNP markers for individual identification of moose (Alces alces), an economically and ecologically important top-herbivore in boreal regions. Reduced representation libraries constructed from 34 moose were high-throughput de novo sequenced, generating nearly 50 million read pairs. About 50 000 stacks of aligned reads containing one or more SNPs were discovered with the Stacks pipeline. Several quality criteria were applied on the candidate SNPs to find markers informative on the individual level and well representative for the population. An empirical validation by genotyping of sequenced individuals and additional moose, resulted in the selection of a final panel of 86 high quality autosomal SNPs. Additionally, five sex-specific SNPs and five SNPs for sympatric species diagnostics are included in the panel. The genotyping error rate was 0.002 for the total panel and probability of identities were low enough to separate individuals with high confidence. Moreover, the autosomal SNPs were highly informative also for population level analyses. The potential applications of this SNP panel are thus many including investigations of population size, sex ratios, relatedness, reproductive success and population structure. Ideally, SNP-based studies could improve today's population monitoring and increase our knowledge about moose population dynamics.

Widespread introgression of mountain hare genes into Fennoscandian brown hare populations.

In Fennoscandia, mountain hare (Lepus timidus) and brown hare (Lepus europaeus) hybridize and produce fertile offspring, resulting in gene flow across the species barrier. Analyses of maternally inherited mitochondrial DNA (mtDNA) show that introgression occur frequently, but unavailability of appropriate nuclear DNA markers has made it difficult to evaluate the scale- and significance for the species. The extent of introgression has become important as the brown hare is continuously expanding its range northward, at the apparent expense of the mountain hare, raising concerns about possible competition. We report here, based on analysis of 6833 SNP markers, that the introgression is highly asymmetrical in the direction of gene flow from mountain hare to brown hare, and that the levels of nuclear gene introgression are independent of mtDNA introgression. While it is possible that brown hares obtain locally adapted alleles from the resident mountain hares, the low levels of mountain hare alleles among allopatric brown hares suggest that hybridization is driven by stochastic processes. Interspecific geneflow with the brown hare is unlikely to have major impacts on mountain hare in Fennoscandia, but direct competition may.

Estimating population size using single-nucleotide polymorphism-based pedigree data.

Reliable population estimates are an important aspect of sustainable wildlife management and conservation but can be difficult to obtain for rare and elusive species. Here, we test a new census method based on pedigree reconstruction recently developed by Creel and Rosenblatt (2013). Using a panel of 96 single-nucleotide polymorphisms (SNPs), we genotyped fecal samples from two Swedish brown bear populations for pedigree reconstruction. Based on 433 genotypes from central Sweden (CS) and 265 from northern Sweden (NS), the population estimates (N = 630 for CS, N = 408 for NS) fell within the 95% CI of the official estimates. The precision and accuracy improved with increasing sampling intensity. Like genetic capture-mark-recapture methods, this method can be applied to data from a single sampling session. Pedigree reconstruction combined with noninvasive genetic sampling may thus augment population estimates, particularly for rare and elusive species for which sampling may be challenging.

Phantoms of the forest: legacy risk effects of a regionally extinct large carnivore.

The increased abundance of large carnivores in Europe is a conservation success, but the impact on the behavior and population dynamics of prey species is generally unknown. In Europe, the recolonization of large carnivores often occurs in areas where humans have greatly modified the landscape through forestry or agriculture. Currently, we poorly understand the effects of recolonizing large carnivores on extant prey species in anthropogenic landscapes. Here, we investigated if ungulate prey species showed innate responses to the scent of a regionally exterminated but native large carnivore, and whether the responses were affected by human-induced habitat openness. We experimentally introduced brown bear Ursus arctos scent to artificial feeding sites and used camera traps to document the responses of three sympatric ungulate species. In addition to controls without scent, reindeer scent Rangifer tarandus was used as a noncarnivore, novel control scent. Fallow deer Dama dama strongly avoided areas with bear scent. In the presence of bear scent, all ungulate species generally used open sites more than closed sites, whereas the opposite was observed at sites with reindeer scent or without scent. The opening of forest habitat by human practices, such as forestry and agriculture, creates a larger gradient in habitat openness than available in relatively unaffected closed forest systems, which may create opportunities for prey to alter their habitat selection and reduce predation risk in human-modified systems that do not exist in more natural forest systems. Increased knowledge about antipredator responses in areas subjected to anthropogenic change is important because these responses may affect prey population dynamics, lower trophic levels, and attitudes toward large carnivores. These aspects may be of particular relevance in the light of the increasing wildlife populations across much of Europe.

Using eDNA to experimentally test ungulate browsing preferences.

Large herbivores may affect ecosystem processes and states, but such effects can be difficult to quantify, especially within multispecies assemblages. To better understand such processes and improve our predictive ability of systems undergoing change, herbivore diets can be studied using controlled feeding trials (or cafeteria tests). With some wildlife, such as large herbivores, it is impractical to empirically verify these findings, because it requires visually observing animals in forested environments, which can disturb them from their natural behaviors. Yet, in field-based cafeteria trials it is nearly impossible to differentiate selection between herbivore species that forage on similar plants and make very similar bite marks. However, during browsing ungulates leave saliva residue which includes some buccal cells and DNA that can be extracted for species identification. Here we used a newly developed eDNA-based method (biteDNA) to test the browsing preferences of four sympatric ungulate species in the wild. Overall, food preferences varied between species, but all species strongly preferred deciduous over coniferous species. Our method allows the study of plant-animal interactions in multispecies assemblages at very fine detail.

Single nucleotide polymorphism-based dispersal estimates using noninvasive sampling.

Quantifying dispersal within wild populations is an important but challenging task. Here we present a method to estimate contemporary, individual-based dispersal distance from noninvasively collected samples using a specialized panel of 96 SNPs (single nucleotide polymorphisms). One main issue in conducting dispersal studies is the requirement for a high sampling resolution at a geographic scale appropriate for capturing the majority of dispersal events. In this study, fecal samples of brown bear (Ursus arctos) were collected by volunteer citizens, resulting in a high sampling resolution spanning over 45,000 km(2) in Gävleborg and Dalarna counties in Sweden. SNP genotypes were obtained for unique individuals sampled (n = 433) and subsequently used to reconstruct pedigrees. A Mantel test for isolation by distance suggests that the sampling scale was appropriate for females but not for males, which are known to disperse long distances. Euclidean distance was estimated between mother and offspring pairs identified through the reconstructed pedigrees. The mean dispersal distance was 12.9 km (SE 3.2) and 33.8 km (SE 6.8) for females and males, respectively. These results were significantly different (Wilcoxon's rank-sum test: P-value = 0.02) and are in agreement with the previously identified pattern of male-biased dispersal. Our results illustrate the potential of using a combination of noninvasively collected samples at high resolution and specialized SNPs for pedigree-based dispersal models.

DNA left on browsed twigs uncovers bite-scale resource use patterns in European ungulates.

Fine-scale resource use by large herbivores is often difficult to quantify directly. This is particularly true for browsing ungulates due to the challenges in observing shy subjects in forested environments of low visibility. As a consequence we know relatively little about resource use by diverse browsing ungulates. When browsing, ungulates leave behind saliva on the browsed twig that includes their DNA, which can be used to identify the species that was responsible for browsing the twig. We used this method, which we term "biteDNA", to study bite-scale browsing patterns in a temperate ungulate community. This approach provides a level of detail in browsing patterns across species that was previously very hard to attain. We found that all deer species largely overlapped in terms of the tree species they used. Moose browsed larger diameters than red deer and roe deer, but these latter two species did not differ. Moose browsed at higher heights than red deer, and red deer higher than roe deer. Although the deer species differed in mean browsing height, species were comparable in terms of their minimum browsing height of ~20 cm. This means that height and diameter ranges of the smaller species were found to be completely inside the ranges of the larger species. Hence, while moose may access exclusive food resources in terms of browse height and diameter, red and roe deer cannot.

De novo SNP discovery in the Scandinavian brown bear (Ursus arctos).

Information about relatedness between individuals in wild populations is advantageous when studying evolutionary, behavioural and ecological processes. Genomic data can be used to determine relatedness between individuals either when no prior knowledge exists or to confirm suspected relatedness. Here we present a set of 96 SNPs suitable for inferring relatedness for brown bears (Ursus arctos) within Scandinavia. We sequenced reduced representation libraries from nine individuals throughout the geographic range. With consensus reads containing putative SNPs, we applied strict filtering criteria with the aim of finding only high-quality, highly-informative SNPs. We tested 150 putative SNPs of which 96% were validated on a panel of 68 individuals. Ninety-six of the validated SNPs with the highest minor allele frequency were selected. The final SNP panel includes four mitochondrial markers, two monomorphic Y-chromosome sex-determination markers, three X-chromosome SNPs and 87 autosomal SNPs. From our validation sample panel, we identified two previously known parent-offspring dyads with reasonable accuracy. This panel of SNPs is a promising tool for inferring relatedness in the brown bear population in Scandinavia.

Browsed twig environmental DNA: diagnostic PCR to identify ungulate species.

Ungulate browsing can have a strong effect on ecological processes by affecting plant community structure and composition, with cascading effects on nutrient cycling and animal communities. However, in the absence of direct observations of foraging, species-specific foraging behaviours are difficult to quantify. We therefore know relatively little about foraging competition and species-specific browsing patterns in systems with several browsers. However, during browsing, a small amount of saliva containing buccal cells is deposited at the bite site, providing a source of environmental DNA (eDNA) that can be used for species identification. Here, we describe extraction and PCR protocols for a browser species diagnostic kit. Species-specific primers for mitochondrial DNA were optimized and validated using twigs browsed by captive animals. A time series showed that about 50% of the samples will amplify up to 12 weeks after the browsing event and that some samples amplify up to 24 weeks after browsing (12.5%). Applied to samples of natural browsing from an area where moose (Alces alces), roe deer (Capreolus capreolus), fallow deer (Cervus dama) and red deer (Cervus elaphus) are sympatric, amplification success reached 75%. This method promises to greatly improve our understanding of multispecies browsing systems without the need for direct observations.

Adaptive suppression of subordinate reproduction in cooperative mammals.

Attempts to account for observed variation in the degree of reproductive skew among cooperative breeders have usually assumed that subordinate breeding has fitness costs to dominant females. They argue that dominant females concede reproductive opportunities to subordinates to retain them in the group or to dissuade them from challenging for the dominant position or that subordinate females breed where dominants are incapable of controlling them. However, an alternative possibility is that suppressing subordinate reproduction has substantive costs to the fitness of dominant females and that variation in these costs generates differences in the net benefits of suppression to dominants which are responsible for variation in the frequency of subordinate breeding that is not a consequence of either reproductive concessions or limitations in dominant control. Here, we show that, in wild Kalahari meerkats (Suricata suricatta), the frequency with which dominants evict subordinates or kill their pups varies with the costs and benefits to dominants of suppressing subordinate breeding, including the dominants' reproductive status, the size of their group, and the relatedness of subordinates. We review evidence from other studies that the suppression of reproduction by subordinates varies with the likely costs of subordinate breeding to dominants.

Factors affecting the reproductive success of dominant male meerkats.

Identifying traits that affect the reproductive success of individuals is fundamental for our understanding of evolutionary processes. In cooperative breeders, a dominant male typically restricts mating access to the dominant female for extended periods, resulting in pronounced variation in reproductive success among males. This may result in strong selection for traits that increase the likelihood of dominance acquisition, dominance retention and reproductive rates while dominant. However, despite considerable research on reproductive skew, few studies have explored the factors that influence these three processes among males in cooperative species. Here we use genetic, behavioural and demographic data to investigate the factors affecting reproductive success in dominant male meerkats (Suricata suricatta). Our data show that dominant males sire the majority of all offspring surviving to 1 year. A male's likelihood of becoming dominant is strongly influenced by age, but not by weight. Tenure length and reproductive rate, both important components of dominant male reproductive success, are largely affected by group size and composition, rather than individual traits. Dominant males in large groups have longer tenures, but after this effect is controlled, male tenure length also correlates negatively to the number of adult females in the group. Male reproductive rate also declines as the number of intra- and extra-group competitors increases. As the time spent in the dominant position and reproductive rate while dominant explain > 80% of the total variance in reproductive success, group composition thus has major implications for male reproductive success.